Distinguishing among dark matter annihilation channels with neutrino telescopes 
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We investigate the prospects for distinguishing dark matter annihilation channels using the neu- 
trino flux from gravitationally captured dark matter particles annihilating inside the sun. We show 
that, even with experimental error in energy reconstruction taken into account, the spectrum of con- 
tained muon tracks may be used to discriminate neutrino final states from the gauge boson/charged 
lepton final states and to determine their corresponding branching ratios. We also discuss the effect 
of v T regeneration inside the sun as a novel method to distinguish the flavor of final state neutrinos. 
This effect as evidenced in the muon spectrum becomes important for dark matter masses above 
300 GeV. Distinguishing primary neutrinos and their flavor may be achieved using multi-year data 
from a detector with the same capability and effective volume as the IceCube/DeepCore array. 
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I. INTRODUCTION 

Many lines of evidence support the existence of dark 
matter (DM) in the universe, but its identity remains a 
major problem at the interface of particle physics and 
cosmology. The weakly interacting massive particles 
(WIMPs) are promising DM candidates [l| • They can ex- 
plain the DM relic abundance, as precisely measured _b 
cosmic microwave background (CMB) experiments 
via thermal freeze-out of annihilation in the early uni- 
verse. Major experimental efforts seek to detect DM 
particles via direct, indirect and collider searches. These 
experiments provide complementary information on the 
properties of DM such as its mass, elastic scattering cross 
section with nucleons, annihilation cross section, and an- 
nihilation channels. 

Indirect searches investigate annihilation of DM to var- 
ious final states (neutrinos, photons, charged particles) 
through astrophysical observations. Neutrinos provide 
an especially interesting probe because they are least af- 
fected on their way from the production point to the 
detection point. As a result, neutrino telescopes like Ice- 
Cube/DeepCore (IC/DC) can trace a neutrino signal di- 
rectly back to the source. DM particles gravitationally 
captured inside the sun annihilate and produce such a 
signal. When equilibrium between DM capture by the 
sun and DM annihilation inside the sun is established, 
the flux of produced neutrinos depends on the DM mass, 
scattering cross section off nucleons and the annihilation 
spectrum. The DM mass and scattering cross section can 
be independently determined from other experiments. In 
particular, the Large Hadron Collider (LHC) is on the 
verge of discovering new physics, which will enable us to 
measure the mass of the DM particle. Therefore, using 
the LHC measurements and the IC/DC results in tan- 
dem, we hope to identify the annihilation channels and 
their corresponding branching ratios. 

In this paper, we study the prospects for determining 
DM annihilation final states with IC/DC. This model- 
independent study has two goals: 1) distinguishing sce- 
narios in which DM annihilates into neutrinos from sce- 
narios where DM annihilates into gauge bosons and 



charged leptons, and 2) extracting information about fla- 
vor composition of neutrino final states in the former 
case. 

Gauge boson and charged fermion final states are the 
dominant annihilation channels for the popular and ex- 
tensively studied neutralino DM in the minimal super- 
symmetric standard model (MSSM). On the other hand, 
annihilation to neutrinos can arise as the dominant chan- 
nel in extensions of the standard model (SM) in which 
DM is related to the neutrino sector. Discriminating neu- 
trino final states from gauge boson and tau final states 
will therefore allow us to discern between the two large 
classes of DM models. In the case that DM mainly anni- 
hilates into neutrinos, it will be important to also know 
the flavor composition of the final state neutrinos. This 
will shed further light on the specific aspects of the model 
that connect DM to neutrinos. 

We show that the contained muon tracks at neutrino 
telescopes like IC/DC may be used for both purposes. 
The spectrum of contained muons can be used to deter- 
mine the branching ratios of neutrino final states versus 
gauge boson and tau final states. We demonstrate this 
for several points in the parameter space within the reach 
of the one-year sensitivity limits of the IC/DC array, with 
the background from atmospheric neutrinos and the ex- 
perimental error in energy reconstruction of the muons 
taken into account. In addition, the v T regeneration due 
to charged current interactions inside the sun may be 
used to extract information about the flavor of final state 
neutrinos for DM masses above 300 GeV. In particular, 
we see that the v T final state can be discriminated from 
the i/fj, and v e final states at a significant level for DM 
masses as heavy as 800 GeV. Distinguishing the neutrino 
final states and discriminating the flavor of neutrinos may 
be achieved with multi-year data from IC/DC. 

This paper is organized as follows. The DM signal in 
neutrino telescopes is explained in section II. In section 
III, we discuss the motivation for models where DM an- 
nihilates mainly into neutrinos. Section IV includes the 
analysis distinguishing the neutrino final state from the 
gauge boson and tau final states. We investigate discrim- 
inating the neutrino flavor in the former case in section 
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V. Finally, we conclude the paper in Section VI. 

II. NEUTRINO TELESCOPES AS DARK 
MATTER DETECTORS 

A. The Sun as a Source of Dark Matter Neutrinos 

DM annihilation in the sun produces a neutrino flux 
that can be probed by neutrino telescopes like IC/DC. 
This neutrino flux is modeled by calculating the number 
of gravitationally captured DM particles in the sun and 
then considering the propagation and detection of the 
neutrinos produced in DM annihilation. The number of 
captured DM particles as a function of time is governed 
by a differential equation that balances the capture of 
particles from elastic scattering on nucleons with the an- 
nihilation rate. The total rate of annihilation in the sun 
is given by 0] (also see (ij) 

r * -§*-*(£)■ « 

where C is the capture rate of DM particles. The number 
of particles captured will saturate at « C/2 as long 
as the length of time for t he p rocess has exceeded the 
equilibration time, T eq = {\f C A)^ 1 . 

With a nominal freeze-out annihilation rate of 3x 10 -26 
cm 3 /s, DM readily achieves equilibrium within the life- 
time of the solar system, 4.5 Gyr, for spin-independent 
elastic scattering cross sections compatible with the lim- 
its set by the XENON100 experiment: cr si = few x 
(10" 9 - 1(T 8 ) pb for the mass range 100 GeV - 1 
TeV The spin-dependent cross section needed to 

reach equilibrium is larger by a factor of about 300 Q, 
which is well below the current experimental limit of 
ctsd < 2.5 x (10" 4 - 10~ 3 ) pb for the 100 GeV - 1 TeV 
mass range Q. 

The difference between the equilibrium constraints on 
<7si and osd comes from the fact that trsi cx M 4 while 
ctsd <x J {J + l)Af 2 ; M and J are the mass and spin of 
target nuclei inside the sun respectively. Heavy elements 
such as iron account for only a thousandth of the sun's 
mass and suffer a two order of magnitude form-factor 
suppression in capturing WIMP-scale DM masses [3j . On 
the other hand, hydrogen dominates the sun's composi- 
tion and does not suffer a form factor suppression. Thus, 
the ratio of the spin-independent contribution from iron 
to the spin-dependent contribution from hydrogen is on 
the order of (56) 4 /10 5 ~ 0(100). Since equilibrium is 
easily achieved in the sun, the neutrino signal will de- 
pend solely on C, or equivalcntly ctsi and/or usd- 1 

Once neutrinos are produced from DM annihilation in 
the sun, their spectra undergo a number of changes be- 



1 The neutrino signal from DM annihilation in the earth is negli- 
gible assuming a standard DM velocity distribution in the halo. 



fore detection. Scattering via neutral current interactions 
results in neutrino energy loss, while charged current in- 
teractions result in absorption of neutrinos. These effects 
are proportional to the energy of the neutrino and will 
therefore have the greatest effect on the high energy part 
of the spectrum. Neutrinos are also affected by oscilla- 
tions on the way to the earth. Since the oscillation length 
is proportional to energy, L osc oc E v , the effect is more 
important for the low energy part of the spectrum. 



B. Neutrino Background, Energy Reconstruction 
and Thresholds 

Neutrino telescopes access the neutrino signal from 
DM annihilation by recording Cerenkov light from rel- 
ativistic charged particles in their volume. Muon neu- 
trinos produce muons via charged current interactions in 
the detector. Cosmic ray showers create a muon back- 
ground that can be controlled by selecting for upward- 
going events since muons are stopped in the earth. This 
limits observation of a DM signal from the sun to half the 
year, when the sun is below the horizon. 2 With through- 
going muons eliminated as a background, the most sig- 
nificant contribution to the remaining background comes 
from atmospheric neutrinos. The spectrum of the atmo- 
spheric neutrino background from cosmic rays is under- 
stood theoretically to within 20% and is measured 
above 100 GeV to within 10% but individual at- 
mospheric neutrino events cannot be distinguished from 
neutrinos from DM annihilation. 

The DM neutrino signal can be further enhanced over 
the background by allowing for an angular cut in the di- 
rection of the sun. The cut will be made on the muon 
track, and not the incoming neutrino. A smaller an- 
gle between the track and incident neutrino occurs for 
muons with an energy closer to that of the neutrino. 
For an elastic collision, the dependence of the angle 
on the muon energy angle is approximately given by 

Out, ^ 5 - 7 ° ( 10 ° GeV/^) 1 / 2 (for example, see @). That 
is, the muon events that deflect little from the incom- 
ing neutrino path are the highest energy events. While 
the smallest possible angular cut is desirable to eliminate 
background, more accommodating cuts on the angle be- 
tween the track and the sun provide information about 
lower energy events. 

Energy reconstruction of events can be approximated 
in two regimes. Above 1 TeV, the light generated 



2 In addition, a portion of the detector may be used as a veto to 
observe contained muon events with a conversion vertex of !/„ 
to n inside the instrumented volume, as in the case of IC/DC. 
The veto procedure virtually eliminates the contribution to the 
background from through-going atmospheric muons by selecting 
for contained vertices. This increases the potential observation 
time to the full year when the sun is both above and below the 
horizon. 
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by an event is proportional to the muon energy since 
Bremsstrahlung radiation, nuclear interactions and pair 
production create charged particles each of which con- 
tribute to the Cerenkov radiation. Energy reconstruction 
in this regime will be logarithmic in energy with an error 
of log 10 <je ~ 0.3 PH. Meanwhile, below - 300 GeV the 
length of the track is proportional to the energy of the 
incident particle since the majority of energy loss in this 
range is governed by ionization. In this case reconstruc- 
tion will likely be more accurate for contained tracks, 
and reconstruction algorithms are currently in develop- 
ment Events between these regimes call for a more 
complicated reconstruction algorithm. 

Whereas the IC design is focused on event energies 
above a TeV, IC /DC achieves an energy threshold as low 
as 10 GeV. IC/DC consists of eight more densely instru- 
mented strings with high quantum-efficiency digital op- 
tical modules (DOMs). Surrounding IC/DC strings veto 
through-going muons so IC/DC records muon- neutrino 
vertex events inside its volume. IC/DC increases the IC 
effective volume at energies below 65 GeV and accounts 
for the majority of events recorded below 100 GeV. Fur- 
ther infills of the IC/DC array, such as PINGU, could 
extend the energy threshold to a few GeV and further 
increase the effective volume by a factor of two at 10 
GeV [3. 



III. PRIMARY NEUTRINOS FROM DARK 
MATTER ANNIHILATION: THEORETICAL 
MOTIVATION 

DM particles can in principle annihilate into any of 
the SM particles. The annihilation rate to a final state is 
given by a ann v = a + bv 2 , where v is the relative velocity 
of annihilating particles. The two terms in this expres- 
sion represent the S-wave and P-wave contributions re- 
spectively. For DM annihilation inside the sun the ther- 
malizcd velocity is very low (v ~ 10~ 4 for a 100 GeV 
DM particle). Therefore channels that proceed through 
the S-wave dominate annihilation (unless a is extremely 
small) . 

For the popular and extensively studied neutralino DM 
in the MSSM, annihilation is mainly into gauge boson 
and charged fermion final states. For example, in the 
stau-neutralino coannihilation region of minimal super- 
gravity (mSUGRA) models [13j . where DM is mostly 
Bino, annihilation in the S-wave of DM particles is typ- 
ically dominated by the tau final states. Taus produce 
relatively soft neutrinos via three-body decay. In the fo- 
cus point region of mSUGRA models [lij ]. where DM 
particles have a large Higgsino fraction, S'-wave annihi- 
lation is predominantly into W bosons (and t-quarks if 
kinematically allowed). W^'s produce neutrinos with a 
harder spectrum via two-body decay. As a result, anni- 
hilation of neutralino DM typically yields secondary neu- 
trinos from W and tau final states. The neutrino signal 
that is produced from neutralino annihilaiton inside the 



sun has been studied in various cases [15J . Because of the 
existence of charged particles in the final state, annihila- 
tion of neutralino DM also results in a gamma-ray signal. 
This implies that one can also obtain constraints by using 
data from the Fermi gamma-ray space telescope [l6| ■ 

Neutralino annihilation to primary neutrinos happens 
via gauge interactions where neutrinos are produced 
through the Zvv vertex. It therefore produces a left- 
handed (LH) neutrino and a right-handed (RH) antineu- 
trino. Considering that the neutralino is a Majorana 
fermion, the S'-wave contribution to such a final state 
is extremely small due to the tiny masses of neutrinos. 
We note that the P-wave contribution to neutrino final 
states is also very small because of velocity suppression. 
In consequence, production of primary neutrinos with a 
hard spectrum from neutralino annihilation is highly sup- 
pressed. 

However, it is possible to enhance DM annihilation 
into primary neutrinos by going beyond neutralino DM. 
If DM annihilation produces a vv pair, instead of a vv 
pair, it can proceed in the S-wave without any mass sup- 
pression. A detailed analysis of settings in which DM 
annihilation into neutrinos is enhanced has been given 
in [l7| . This can happen when DM is related to the neu- 
trino sector (for some specific models, see [HI)- 

One particular model that can produce hard 
monochromatic neutrinos is the B — L extension of the 
MSSM (B and L are baryon and lepton number respec- 
tively). A minimal extension of the SM gauge group, 
motivated by the nonzero neutrino masses, includes a 
gauged U{1)b-l gauge symmetry [19(. Anomaly cancel- 
lation then implies the existence of three RH neutrinos 
(N), the lightest superpartner of which (TV) makes a vi- 
able DM candidate |20j ■ Thermal freeze-out of the sneu- 
trinos can yield the correct dark matter abundance if the 
B — L symmetry is broken around the TeV scale [2lj . 
The sneutrinos can annihilate in the S-wave to RH neu- 
trinos, which in turn decay to LH neutrinos and the SM 
Higgs. The energy of the resulting neutrinos will be close 
to monochromatic as long as the difference between the 
masses of the RH sneutrino and RH neutrino is much less 
than the RH sneutrino mass. The prospects for DM de- 
tection for this model using IC/DC is considered in 22[. 3 

The neutrino signal is the main channel of indirect 
searches for models with DM annihilation to primary 
neutrinos. These models result in a highly suppressed 



3 This model has other interesting cosmological and phenomeno- 
logical implications. The Tevatron \2i and LEP |25H limits on 
the Z' mass bound the RH sncutrino-nucleon scattering cross 
section to be a < 7 X 10~ 9 pb, which is just below the bound 
from the XENON experiment [5(. Significant DM annihilation to 
taus is also possible 20] (also see [2611 ). which could provide an 
explanation for the positron excess in the cosmic rays reported 
by PAMELA |27|1 . In addition, this model can accommodate in- 
flation 28] , provide a unified picture of inflation and dark mat- 
ter [2ll . and give rise to interesting predictions for neutrinoless 
double beta decay experiments |29( . 
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gamma-ray signal, which escapes the bounds set by 
Fermi [l6]. It is worth pointing out that neutrinos from 
DM annihilation in the galactic center provide comple- 
mentary information to those from annihilation inside 
the sun: the former constrains the DM annihilation cross 
section (23|, while the latter bounds the DM scattering 
cross section off nucleons |6J . The signature of DM an- 
nihilation to neutrinos at the galactic center should be 
distinctive with a hard cutoff (for example, see (30jV 



IV. DISTINGUISHING NEUTRINO FINAL 
STATE FROM W AND TAU FINAL STATES 

MSSM neutralino annihilation predominantly pro- 
duces gauge bosons and charged fermions, each of which 
in turn produce secondary neutrinos via two-body or 
three-body decays. In this work, we adopt a model- 
independent approach and consider interesting neutrino 
spectra from DM annihilation to prompt u'a, Ws, and 
r's. 

For charged fermions, we investigate only annihilation 
to r's for the following reasons, e's and fi's are stopped 
immediately in the sun, so [i decay will occur at low en- 
ergy and result in a very soft spectrum. When DM an- 
nihilation produces quarks, all quarks except the i-quark 
hadronize before their subsequent decay. The i-quark de- 
cays to a W boson and a 6-quark before hadronization, 
and the subsequent W decay produces a neutrino spec- 
trum comparable to that from the W final state. For the 
remaining hadronizing quarks, the lighter the quark the 
longer it will take to decay, and hence the lower its energy 
will be. The 6-quark produces neutrinos via three-body 
decay, which has a softer spectrum than that from t. All 
other quarks decay effectively at rest and have unmea- 
surable spectra at low energies that cannot compete with 
the atmospheric background [3]. 

With regard to gauge boson final states, W and Z spec- 
tra are relatively similar; we consider only W final states 
in our analysis. Finally, the Higgs boson final states also 
result in a soft neutrino spectrum that is negligible com- 
pared to that from W and r final states. The Higgs 
mainly decays to 6-quarks, thus yielding a neutrino spec- 
trum similar to that from the 6-quark final state. 

In the following, we therefore focus on secondary neu- 
trinos from DM annihilation to W and r final states vs 
primary neutrinos and their corresponding spectra. Pre- 
vious investigations have shown that measuring the spec- 
trum can allow the reconstruction of the DM mass and 
its annihilation branching ratios H (a lso see [3l| for re- 
constructing DM properties, and [32| for the signature 
of primary neutrinos). We perform an analysis for re- 
constructing the DM elastic scattering cross section and 
its annihilation channels, including experimental error in 
energy reconstruction and assuming that its mass can 
be determined from other experiments (notably by the 
LHC). 



A. Neutrino & Contained Muon Spectra 

In this analysis, we consider the spectrum of contained 
muon tracks so that the vertex at which is converted 
to a fi is within the detector volume. Hadronic and elec- 
tromagnetic cascade events in the ice also carry useful 
information about low energy neutrinos and may show 
tau neutrino appearance from oscillations [331 ]. IC/DC is 
uniquely suited to measure these low energy events and 
recently confirmed the observation of neutrino-induced 
cascades [121 ] . For our purposes, since cascades are local- 
ized in the ice and thus suffer from an angular resolution 
of about 60 degrees at 100 GeV, too much background is 
admitted to perform a meaningful DM search with cas- 
cades [aj . 4 

We use DarkSUSY to calculate the spectrum for any 
given final state (36|. To the extent that the energy of 
the event may be well reconstructed in this range, we 
also assume that the muon track itself is fully contained, 
ending inside the detector. We do not account for any 
loss of events or poor energy reconstruction at the edges 
of the fiducial volume. 

As an example, Fig.Q]depicts the theoretical prediction 
for neutrino and muon spectra that result from the an- 
nihilation of 100 GeV DM particles in the sun to prompt 
Vp's, W's and r's. The flux is given in events per km 3 
effective volume of the detector at earth per year. Only a 
half-year of data is taken for each calendar year, since we 
assume the neutrino telescope requires the sun to be be- 
low the horizon to observe the DM neutrino signal above 
the atmospheric muon background. For a DM mass of 
100 GeV the muons will be fully contained and accessi- 
ble in IC. The qualitative features of the neutrino and 
muon spectra (such as peaks and kinematic cutoffs) are 
the same for higher DM masses. 

The prompt neutrinos, Fig. Hal result in a monochro- 
matic peak at the DM mass in the neutrino spectrum. 
The contained muon spectrum is plotted along with the 
neutrino spectrum to emphasize the effect of the charged 
current conversion of neutrinos to muons. The conversion 
demonstrates the weak scale reduction in events relative 
to the incident neutrinos and is proportional to energy, 
so the peak from the neutrino annihilation channel is 
well preserved. The spectra are largely insensitive to the 
choice of neutrino oscillation parameters and annihilating 
neutrino flavor. Here we depict 813 = 10°, normal mass 
hierarchy and annihilation to f M , but scenarios with in- 
verted mass hierarchy, < 813 < 10° and other neutrino 
flavors in the final state only slightly change the height 
of the peak. 5 



4 Muon tracks and cascades together can be used to gain informa- 
tion about the flavor ratio of neutrinos at the detector, which 
can be used as diagnostic of DM annihilation inside the sun [35|l . 

5 Recent results from T2K [13 and MINOS ^ experiments sug- 
gest nonzero 6*13 at 2.5ct and 89% confidence respectively. The 
results are consistent with 613 = 10°. 
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(b) Annihilation to W final state 
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(c) Annihilation to r final state 
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This is the reason for the relatively sharp kinematic edges 
in the neutrino spectrum in Fig. Ilbl We note that W's 
also decay to 6-quarks, which in turn produce tertiary 
neutrinos. This results in additional contributions to the 
neutrino spectrum below the lower kinematic cutoff. The 
muon spectrum has a peak at energies well below the DM 
mass resulting from charged current interactions whose 
cross section is proportional to the neutrino energy. 

An important point to note is that W's with transverse 
polarization can dominate the annihilation final state as 
happens, for example, in the case of neutralino DM. The 
main contribution to hard neutrinos (with an energy close 
to the upper kinematic cut off mentioned above) comes 
from transverse W's in this case, which will become to- 
tally dominant when Mdm ^> My/. In these cases, it 
is therefore important to carry out a spin-dependent cal- 
culation that accounts for spin correlations of the final 
state particles and the helicity- dependence of their de- 
cays [30| . We also note that the upper kinematic cut 
off itself gets closer to Mdm, and thus creates a harder 
neutrino spectrum, when Mum ^ Mw- 

The t's produce neutrinos via three-body decays, 
which results in a much softer neutrino spectrum that 
rises toward lower energies as seen in Fig. [Tc] Thus the 
peak of the resulting muon spectrum is located at a lower 
energy than that in Fig. Ilbl However, there is no kine- 
matic cutoff in the neutrino spectrum in this case, and 
hence the muon spectrum extends smoothly to the DM 
mass. 

The important point is that there is a distinct sep- 
aration in energy between the contained muon peak of 
prompt neutrinos at the DM mass and the W and r peaks 
at lower energies. This indicates that scenarios with DM 
annihilation to primary neutrinos can be distinguished 
from those with DM annihilating to W bosons or t's. 
Moreover, were a model to contain both a neutrino final 
state and a W or r final state, the neutrino signal may 
be used to determine the corresponding branching ratios. 



FIG. 1: Energy spectra of contained muon events (solid) and 
incident neutrinos (dashed) at the earth from a 100 GeV DM 
particle for various annihilation channels. Relative to the sun 
1°, 5° angular cuts have been placed on the signal and atmo- 
spheric background muons (thin black lines) . The top plot for 
annihilation to neutrinos is for ctsd = 10 -41 cm 2 , while the 
other plots are for aso = 10~ 40 cm 2 , all of which are below 
the current experimental bounds. Normal mass hierarchy and 
#13 = 10° have been chosen here, but the spectra are largely 
insensitive to the oscillation parameters. 



The W bosons produce neutrinos via two-body decay, 
which are softer than the previous case. For a highly 
boosted W (whose energy is equal to the DM mass), the 
energy of secondary neutrinos is between Mom{^— /?)/2 
and M DM (1 + 0)/2, where /3 = y/1 - (M W /M DM ) 2 - 



B. Reconstruction of Annihilation Channels 

The theoretical predictions shown in the previous sub- 
section are subject to experimental error. Individual 
muon events from a particular annihilation channel can- 
not be tagged as such, nor can individual signal and back- 
ground events be distinguished. A more realistic picture 
is shown in Fig. [2] for a 100 GeV DM particle, where two 
channels along with background events have been added 
together and subjected to a 5° angular cut. The branch- 
ing ratios are 90% to W bosons and 10% to neutrinos 
with ctsd = 10~ 40 cm 2 . This is a factor of a few below 
the current bounds from IC [|| and just at the current 
bounds of Super Kamiokande (SuperK) for annihilation 
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to W's d3. 6 

Angular cuts on the muon tracks relative to the posi- 
tion of the sun lower the relevant atmospheric contained 
muon background relative to the signal. Lower energy 
events are preferentially scattered at higher angles rel- 
ative to the incoming neutrino; therefore, angular cuts 
disproportionately remove low energy events. A 1° cut 
leaves the majority of the high energy muons from pri- 
mary neutrinos, which are clearly separated from a re- 
duced background. Meanwhile, a 5° cut admits more 
lower energy muons that are indicative of the pronounced 
peaks from W and t hnal states. Different angular cuts 
will optimize the signal to background ratio for differ- 
ent annihilation channels. We have found a 5° cut to 
be optimal, where the background does not overwhelm 
the signal in the regime above 40 GeV and a significant 
portion of the W peak remains. The contribution from 
neutrinos is still easily distinguished in the peak cutoff. 




20 40 60 80 100 120 140 

E [GeV] 



FIG. 2: Total contained muon spectrum from 90, 100 and 110 
GeV DM particles (from left to right) with ct S d = 10~ 40 cm 2 
annihilating to W's (90%) and v's (10%) with the atmospheric 
background added. A 5° cut has been placed on the events. 
Smeared spectra for energy reconstruction errors of 10 GeV 
(dashed) and 40 GeV (dotted) are also shown. 

The realistic energy reconstruction for a contained 
muon spectrum at IC/DC for these energies is not 
certain. Energy reconstruction for through-going TeV 
muons is approximately 0"i og e = 0-3 [13] ■ Energy loss 
of the muon is described by 

f-. + W, (2) 

where a quantifies muon loss via ionization, and b quan- 
tifies loss from pair productions, Brehmsstrahlung radi- 
ation and nuclear interactions. For ice at energies above 



For a purely spin-independent cross section, the same result is 
obtained for crgi S3 3 X 10 — 43 cm 2 . However, such a large value 
of (Tgj is already ruled out by direct detection experiments 
For this reason we focus on the spin-dependent cross sections in 
this analysis. 



1 TeV, enough Brehmsstrahlung and pair production oc- 
curs to ensure that the light deposited in the detector is 
proportional to the muon energy. However, for contained 
muon events less than one kilometer in extent, the track 
length is proportional to the muon energy. Reconstruc- 
tion in this regime will be better than in the regime above 
1 TeV but will depend on the geometry of the detector. 

The energy resolution error of the IC/DC detector for 
contained muon events about 100 GeV should be linear 
with energy, in that the track length is proportional to 
the energy of the incoming muon. In Fig. [21 we recreate a 
Gaussian energy reconstruction with an error, or width, 
equal to 10 or 40 GeV. This is equivalent to claiming that 
the track length can be known to within 50 — 200 m, a 
resolution that seems feasible for IC/DC (or telescopes 
withy similar capabilities) since IC/DC DOMs are spaced 
at 7 m vertically and strings are 72 m apart. 

The experimental error in measuring the DM mass at 
colliders should also be taken into account. It has been 
shown that in some supersymmetric scenarios a 10% ac- 
curacy can be reached at the LHC For a 100 GeV 
DM particle this corresponds to an error of 10 GeV, and 
hence we include 90, 100 and 110 GeV masses in our 
analysis. 

It is seen that the smoothing of the spectrum consider- 
ably suppresses the distinct channel features, especially 
for a 40 GeV error in energy reconstruction. Never- 
theless, it is still possible to reconstruct the total spin- 
dependent scattering cross section ctsd and the branching 
ratio for annihilation to primary neutrinos BR^ . We also 
note that the uncertainly in determination of the DM 
mass is subdominant to the experimental error in energy 
reconstruction, which results in minimal changes in the 
smoothed spectra. 

In Fig.[3]we show the 2a confidence contours for 90, 100 
and 110 GeV DM particles annihilating to v and W final 
states, using one year of data and applying energy recon- 
struction errors of 10 GeV and 40 GeV. These contours 
account for Poisson errors, the addition of background 
and a 40 GeV energy cut; we have used the x 2 analysis 
in [4l[ to obtain these results. We again see that varying 
DM mass from 90 GeV to 110 GeV produces minimal 
changes in the contours. 

The plot in Fig. [3] also shows the current SuperK 
bound [12] and the future IC/DC sensitivity bound [23[ 
on (Tsd (dashed and dotted curves respectively). For 
100% DM annihilation to W's, the experiments have de- 
rived bounds on osd by applying appropriate energy and 
angular cuts on the muon spectrum. No limit has been 
placed yet for DM annihilation to primary neutrinos by 
either experiment. In the absence of a thorough analysis, 
we have used a simple criterion to find an approximate 
bound on ergo for the neutrino final states. We compared 
the total number of contained muon events with a 40 GeV 
energy threshold and 5° angular cut for the two final 
states. In the case of primary neutrinos, the flux of /it's is 
about eight times larger than the flux from W's, as ex- 
pected from the harder neutrino spectrum. We assumed 
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FIG. 3: Contours of 2a confidence for reconstructing spin- 
dependent cross section ctsd and branching ratio to prompt 
neutrinos BR„ for DM particle annihilating to i/'s and W'a. 
Contours for 90 (inner), 100 (middle) and 110 GeV (outer) 
DM are shown. One year of data with km 3 effective volume 
is used with energy reconstruction errors of 10 GeV (top) 
or 40 GeV (bottom). The region of interest is between the 
current SuperK upper bound (dashed) and the future 
IC86 sensitivity bound (dotted) HI]. 



the limit on <tsd in this case is also tighter by the same 
factor. A dedicated analysis for the neutrino final state, 
with optimized angular and energy cuts, would result in 
a more precise bound. Indeed, theoretical motivation for 
models with DM annihilation to primary neutrinos war- 
rants such a study by experimental collaborations. In 



the presence of both v and W final states, the limits are 
approximately given by 

a SD (1 + 7 BR„) < (i max . (3) 

The current bound from SuperK [42j and the future 
IC/DC sensitivity bound [23| correspond to <7 max « 
2.5 x 10~ 40 and 4 x 10~ 41 cm 2 respectively. 

It is seen from the shape of the contours that higher 
ctsd and smaller BR„ are difficult to distinguish from 
lower (Jsd and larger BR„. The number of muon events 
at the DM mass mainly comes from the peak of primary 
neutrinos, which is determined by BRj,. This holds af- 
ter smoothing of the spectrum because of the kinematic 
cutoff that appears well below the DM mass for annihila- 
tion to W's, see Fig.llbl This implies that after smearing 
the main contribution to the total spectrum at energies 
around and above the DM mass still comes from primary 
neutrinos. On the other hand, the main contribution to 
lower energy muons comes from secondary neutrinos pro- 
duced by W decay. For small values of <tsd or large val- 
ues of BR„ , this contribution is small and overwhelmed 
by the Poisson error from the atmospheric background. 
Therefore, one can compensate for a change in <tsd by a 
corresponding change in BR„ and obtain spectra that are 
statistically indistinguishable. For larger values of ctsd 
or smaller values of BRj, , the contribution of secondary 
neutrinos becomes significant, which limits simultaneous 
variations in ctsd and BR„ that keep the peak height un- 
changed. As a result, the 2a contours are tighter for re- 
construction points toward the right and bottom of the 
plot in Fig. [31 As expected, an energy reconstruction 
error of 40 GeV does a poorer job because it further sup- 
presses the features discussed above. 

In Fig. [5] we show the 2a confidence contours from 
a similar analysis for annihilation to z/s and r's. The 
contours are narrower in this case, implying that the re- 
construction is easier. The main reason for this is that 
the number of muon events from secondary neutrinos is 
significantly larger, compare Figs. [TBI and [Tel The larger 
signal to background ratio places a stronger constraint 
on the simultaneous variation of BR^ and ctsdj which 
results in an improved reconstruction. 

Some comments are in order. In producing these fig- 
ures, a scan of angular cuts in 1° increments within the 
1° — 10° range was made. While visually a 1° cut may 
be optimal in distinguishing the presence of the neu- 
trino channel, a 5° cut was optimal in reconstructing the 
branching ratios in the presence of the broad, low energy 
W or r spectrum. 

We also note that the muon event rates for both the sig- 
nal and background must be convolved with the effective 
volume of the detector for contained muon events, which 
is a function of energ y. I C /DC maintains significant vol- 
ume above 10 GeV [13] . While the geometric volume 
for IC/DC is about 3% of the volume of IC, the more 
dense spacing of DOMs in IC /DC make it more efficient 
at event detections for energies below 100 GeV [l2j. The 
effective volume in this range increases with energy as 
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FIG. 4: The same as Fig. [3] but for annihilation to v and r 
final states. 



longer muon tracks are more likely to produce detectable 
light. This can make the effect of the peak at the DM 
mass from primary neutrinos even more distinctive above 
the smeared W or r spectrum. The reconstruction can 
therefore be better than shown here after accounting for 
the energy dependence of the effective volume. 

Finally, one year results for a km 3 effective volume, 
with the sun below the horizon, translate to roughly ten 
years of a 0.05 km 3 detector capable of operating for a full 
year. Thus, in the region of IC/DC sensitivity that would 
yield a discovery of prompt neutrinos after one year, 10 
years of data could allow reconstruction of the branching 
ratios for the IC/DC effective volume. A larger detector 



with the same capabilities can significantly improve these 
prospects. 



V. DISTINGUISHING NEUTRINO FLAVORS 

Once the presence of direct DM annihilation to neutri- 
nos is confirmed, it is desirable to also learn the flavor of 
the final state neutrinos. This will provide specific infor- 
mation about models that connect DM to the neutrino 
sector. For example, in the U(1)b-l extension of MSSM 
where the lightest RH sneutrino is the DM candidate |20j , 
flavor composition of final state neutrinos is related to the 
neutrino Yukawa couplings. Therefore, knowledge of the 
neutrino flavors will yield useful information pertaining 
to the underlying neutrino mass model and leptogenesis. 

Here we show how regeneration of v T inside the sun and 
its effect on the muon spectrum may be used as a novel 
method to distinguish the flavor of primary neutrinos in 
the final state. 7 



A. The v T Regeneration Effect 

DM annihilation produces neutrinos in flavor eigen- 
states in the sun. Neutrinos then undergo charged cur- 
rent interactions with matter as they propagate through 
the sun. These interactions convert v e , , i/ r to e, /z, r 
respectively, e's and //'s are stopped immediately due to 
electromagnetic interactions. On the other hand, the r 
decays quickly before losing too much energy because of 
its very short lifetime of 3 x 10~ 13 s This decay pro- 
duces a v T , which has a lower energy than the original 
one. Charged current interactions therefore suppress the 
peak of the neutrino spectrum at the DM mass for all fla- 
vors. 8 However, in the case of the u Tl the regeneration of 
neutrinos via three-body decay populates the spectrum 
at energies well below the DM mass. 

The cross section for charged current interactions is 
proportional to the energy of primary neutrinos produced 
at the center of the sun, which is essentially equal to 
the DM mass. Neutrino absorption becomes significant 
when the absorption length of neutrinos L a bs is roughly 
equivalent to the core size of the sun Rq ~ 70, 000 km. 
Using the charged current neutrino-nucleon cross sec- 
tion [44j and a core density of 150 g/cm 3 , we find that 



7 It has been proposed that seasonal variation of the neutrino sig- 
nal may also be used to extract information about the flavor of 
primary neutrinos |45| . 

8 Neutrinos also have neutral current interactions with matter in- 
side the sun. Scatterings via neutral current interactions result 
in energy loss of the neutrinos and further suppress the peak at 
the DM mass. However, the cross section for neutral current in- 
teractions is a factor of 3 smaller than that for charged current 
interactions, which makes them subdominant. More importantly, 
neutral current scatterings affect all flavors equally. 
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FIG. 5: Contained muon spectra for DM annihilation to v £ 
(dotted), Vp, (solid) and v T (dashed) for 400 (Fig. [5a} and 600 
GeV (Fig. HE} DM masses and a S u = 10" 41 cm 2 . Normal 
mass hierarchy and #13 = 10° are chosen here, but substan- 
tive features of the spectra do not depend on the neutrino 
oscillation scenario. A 5° cut has been placed on the events. 
Corresponding neutrino spectra are shown in thin lines. The 
regeneration effect is evident for the v T final state and results 
in a peak in the muon spectrum at low energies. For a 400 
GeV DM particle, regeneration also affects the final state 
due to efficiency of — v T oscillations inside the sun. 



L a bs ~ 70, 000 km at energies E v ~ 300 GeV. Oscillations 
among different flavors should be taken into account as 
neutrinos travel through the sun. We note that for v e 
the flavor and mass eigenstates are the same deep inside 
the sun. This implies that upon production at the cen- 
ter of the sun v e J s propagate through the core without 
changing their flavor. Thus absorption via charged cur- 
rent interactions start to suppress the flux of v e for DM 
masses above 300 GeV as discussed above. 

On the other hand, and v r are not mass eigenstates 
inside the sun, and hence undergo oscillations. Since 
matter effects are the same for these two flavors, the 
Vft — v r oscillation length inside the sun is set by the at- 
mospheric mass splitting L osc = i? ;y /47rAm^ tm . As long 
as L a b s <; L osc /4, oscillations mix and v T efficiently 
before the absorption becomes important. As a result, 



final states also feel the regeneration effect. Since 
L a bs oc E~ x and L osc oc E v , at sufficiently high ener- 
gies Labs drops below L osc /4. This happens for a DM 
mass of about 500 GeV. Starting at this point, — v r 
oscillations cease to be effective. Hence gets absorbed 
through charged current interactions similar to v e . In 
consequence, only the v r final state retains a significant 
regeneration signature for DM masses above 500 GeV. 

In Fig. [5] we show the muon and neutrino spectra 
for different flavors of primary neutrinos. As we see in 
Fig. [5aJ the channel shows some regeneration effect for 
a 400 GeV particle, which makes it distinguishable from 
the v e channel. For a 600 GeV DM particle, see Fig. I5bl 
the v il — v T oscillations are inefficient. Therefore, only the 
v T channel shows significant regeneration, which makes 
it distinguishable from the and v e channels. The re- 
generation effect results in a peak in the muon spectrum 
at low energies, which becomes more pronounced as the 
DM mass increases. These figures are for the normal hi- 
erarchy of neutrinos and #13 = 10°, but variations in the 
neutrino oscillation scheme do not change the substantive 
features of the spectra. 



B. Signals of Neutrino Flavor 

As we saw, the v T regeneration effect becomes signifi- 
cant at DM masses above 300 GeV. The background from 
atmospheric neutrinos is relatively small at such energies 
due to the power law decrease in cosmic ray background. 
However, the neutrino signal from DM annihilation is 
also kinematically suppressed for heavier DM masses. 
This leaves few events to detect in the interesting range of 
the spectrum where a regeneration peak is evident. Ad- 
ditionally, energy reconstruction for the contained muon 
spectrum above 300 GeV begins to suffer from logarith- 
mic error in energy since the reconstruction also depends 
on the amount of light produced in Brehmsstrahlung ra- 
diation and pair production. Events cannot be fully con- 
tained at these energies, so reconstruction efforts based 
on track length are imprecise. 9 

To mitigate these experimental challenges, we inte- 
grate the contained muon events above a threshold to 
consider if the cumulative effect of the regeneration is 
visible. In Fig. [5] we show the integrated events above 60 
GeV with a 5° angular cut, which readily retains the low 
energy regeneration effect. One can see the separation 
between the v e channel and the v^, v T channels above 
300 GeV, and between the v T channel and the v ei 
channels above 500 GeV (as discussed above). Separation 
between signals from different channels is typically larger 
than the Poisson error of the background and signal to- 



If the energy of the corresponding cascades accompanying the 
muon event can also be captured, energy reconstruction would 
be improved. 
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FIG. 6: Integrated muon spectra from 60 GeV to the DM 
mass with a 5° angular cut as a function of DM mass for 
<7SD — 10 -40 cm 2 . The dotted, solid and dashed lines denote 
the u e , and v T final states respectively. The background is 
shown in black line. The bands of lines for each flavor depict 
different oscillation scenarios (normal and inverted mass hier- 
archy as well as < #13 < 10°). The separation between the 
v T and exceeds the Poisson error of the signal and back- 
ground for DM masses above 300 GeV, and the v T becomes 
distinguishable from above 600 GeV as well. 

gether. The background and signal are shown separately 
since the background may be subtracted from the signal 
by observing away from the sun, off-source, as is done 
in galactic center DM searches [23|]. We also note that 
the oscillation scenario does not significantly affect these 
results. The bands of lines depicting the normal and in- 
verted hierarchies as well as #13 ranging from to 10° for 
each flavor do not overlap at DM masses above 300 GeV. 
These results assume ctsd = 10~ 40 cm 2 , which is compat- 
ible with the current bounds from IC for heavier DM, 
and is within the reach of IC/DC sensitivity limits [23j]. 
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FIG. 7: The same as in Fig. [6] after normalization to account 
for unknown value of the DM elastic scattering cross section. 
The normalized spectra give the same number of events (when 
integrated from 200 GeV to the DM mass for a 1° angular cut) 
as the Vfj, final state for normal mass hierarchy and #13 = 10° . 
The separation between the v T and other flavors still exceeds 
the Poisson error of the signal and background. 

If the elastic scattering cross section of the DM is 



known, then a simple integration as in Fig.[B]will be suf- 
ficient to distinguish the flavor of primary neutrinos by 
using the regeneration effect. However, a confirmation 
of annihilation to primary neutrinos would likely come 
from the presence of a monochromatic peak alone, the 
height of which depends on the scattering cross section. 
In the absence of information about the cross section, 
we then need to normalize the signal to extract informa- 
tion about the neutrino flavor. The height of the peak 
in the muon spectrum is closely related to the height of 
the neutrino peak. However, for DM masses above 300 
GeV, for which the regeneration effect becomes impor- 
tant, muon events at the peak are not fully contained. It 
will be most difficult to reconstruct the energy of these 
events without track-length information. Instead, we in- 
tegrate muon events above 200 GeV, which are essentially 
the through-going muons, using a 1° angular cut on the 
muon spectrum in order to best capture the monochro- 
matic peak. We then use this to normalize the integrated 
muon flux above 60 GeV with a 5° cut to retain the max- 
imum effect of regeneration at low energies. 

The normalization accounts for the lack of knowledge 
in the value of the DM elastic scattering cross section, 
and to a lesser extent the oscillation parameters. After 
normalization, see Fig. [7j the separation between the v T 
channel and the v e , channels remains and is still larger 
than the statistical error of the background and signal 
together. Using the regeneration effect, the IC/DC ef- 
fective volume could reasonably yield a determination of 
the neutrino flavor after 10 full years of operation. 



VI. CONCLUSION 

We have investigated prospects of determining DM an- 
nihilation final states with neutrino telescopes by using 
the spectrum of contained muon tracks from conversion 
of neutrinos that are produced in the annihilation of DM 
particles trapped inside the sun. Our focus was on distin- 
guishing neutrino final states from gauge boson and tau 
final states and on discriminating the flavor of final state 
neutrinos. Gauge boson and tau final states are typi- 
cally the dominant annihilation channels in supersym- 
metric models with neutralino DM, while direct annihi- 
lation into neutrinos can occur in models that connect 
DM to the neutrino sector. The theoretical motivation 
for the latter could provide the grounds for a dedicated 
analysis by the IceCube Collaboration to put stringent 
bounds on annihilation to primary neutrinos, similar to 
what has been done for annihilation to gauge bosons [f| . 

Primary neutrinos from DM annihilation result in a 
distinct peak in the muon spectrum at the DM mass. 
For DM masses below 300 GeV we can expect that the 
peak will be accessible to a detector the size of IC. The 
spectrum is smeared as a result of the experimental error 
in energy reconstruction, but primary neutrinos may be 
distinguished from gauge boson and tau final states after 
this effect is taken into account. We showed that for 



11 



an energy resolution of 10 GeV (as in IC/DC) and by 
making an optimal angular cut on the muons (which we 
found to be 5°), the branching ratios may be determined 
in the parameter space within the reach of the one-year 
sensitivity limits of IC /DC with a km 3 /yr of data. This is 
roughly equivalent to 10 years of data for a detector with 
the same capabilities and effective volume as IC/DC. 

The regeneration of v T inside the sun may be used to 
distinguish the flavor of final state primary neutrinos. 
This effect becomes important for DM masses above 300 
GeV and populates the spectrum with muons whose en- 
ergy is well below the energy of primary neutrinos. For 
DM masses up to about 500 GeV, oscillations mix and 
v T effectively, which implies that regeneration affects fi- 
nal states with and v T similarly. Final states with i/ e 
are therefore distinguishable within this mass range. For 
heavier DM particles, the — v r oscillation becomes in- 
efficient. As a result, v r final states are picked out by the 
regeneration effect for DM masses above 500 GeV. We 
showed that final states with v T stand out at a statisti- 
cally significant level for DM masses as heavy as 800 GeV, 



even after normalizing the muon spectrum to the total 
event count (to account for the unknown DM-nucleon 
elastic scattering cross section). Again, such a distinc- 
tion may be achieved with 10 years of data from IC/DC. 

In summary, using the the IC/DC results in tandem 
with independent measurements of the DM mass (for ex- 
ample, from the LHC), will allow us to identify the anni- 
hilation channels of DM with multi-year data. Improved 
energy resolution and increased effective volume of the 
detector will greatly help in achieving this goal. 
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